Exact solution for the critical state in thin superconductor strips with field dependent 

or anisotropic pinning 
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An exact analytical solution is given for the critical state problem in long thin superconductor 
strips in a perpendicular magnetic field, when the critical current density jc{B) depends on the 
local induction B according to a simple three-parameter model. This model describes both isotropic 
superconductors with this jc{B) dependence, but also superconductors with anisotropic pinning 
described by a dependence jc{d) where 9 is the tilt angle of the flux lines away from the normal to 
the specimen plane. 



PACS numbers; 74.60.-w, 74.60.Ge, 74.60.Jg 



I. INTRODUCTION 

The critical state modelEI for the magnetic behavior of 
superconductors with flux-line pinning has proven very 
usefuB though it originally was applied to the simple 
(demagnetization-free) longitudinal geometry of long su- 
perconductors in parallel magnetic field, ft took over 
30 years until an analytical solution of the critical state 
model was obtained for the more realistic transverse ge- 
ometry of thin supierconductoKS. The solutions were de- 
rived for thin diskg3 and stripgj in a perpendicular mag- 
netic field, extending an earliec, work on superconductor 
strips with transport current J3 and finally for elliptic- 
shaped plateletSrfl Recent detailed numerical work for 
stripsQ and disksO of finite thickness shows how the tran- 
sition from longitudinal to transverse geometry occurs 
with changing aspect ratio of the specimen. 

So far, in the transverse geometry all analytical solu- 
tions of the critical state model were restricted to the 
Bean model of constant critical current density jc = 
const, but in many experiments jc = jc{B) depends on 
the local magnetic induction B. For example, the simple 
Kim vaodemjc{B) = jc(0)/(l + \B\/Bo) was considered 
in many experimental and theoretical papers, see e.g. the 
reviews p|,|lO ,p] and the-partly analytical calculations for 
thin stripsE3 and disksEJ. While numerical computa 
easily allow us to consider any jc{B) dependence,LfQlli 
an exact analytical solution of some model may give 
deeper insight since it yields explicit dependences of the 
resulting quantities on the input parameters. 

In the highly anisotropic high- Tc superconductors the 
flux-line pinning in general depends on the angle 9 be- 
tween the local direction of the magnetic induction B 
and the c axis, which in typical experiments is normal 
to the plane of the sample. For example, this type of 
anisotropy occurs when one takes into account the in- 
trinsic pinning exerted by the CuO planes or the pinning 



by extended defects. It has been shown recentlyEl" 
that for thin superconductors of any shape (with thick- 
ness d much smaller than the lateral extension L but 
larger than the magnetic penetration depth A) any such 
out-of-plane-anisotropy of pinning is equivalent to an in- 
duction dependence of the critical sheet current Jc{B) 
(the sheet current is defined as the current density inte- 
grated over the film thickness). Thus, the description of 
the two-dimensional critical state, e.g., in an anisotropic 
strip can be reduced to the analysis of a one-dimensional 
problem with some Jc{B). In this case the characteris- 
tic scale Bq over which Jc{B) changes is of the order of 
1^0 jcd. 

In this paper we present an analytical solution for the 
critical state in thin superconductor strips in perpendic- 
ular field with field dependent critical current density 
jc{B) or, equivalently, with anisotropic pinning described 
by a jc{9)- A three-parameter model jc{B) consisting of 
two straight lines, an inclined line at small B and a hori- 
zontal line at larger B, is considered. This rather general 
model is equivalent to a piecewise constant angular de- 
pendence jc{d) = jci ior < 9 < 00 and jc{9) = jc2 for 
9o < 9 < Tr/2 where jci, jc2, and are the parameters 
of the model. We shall show below that the steepness of 
the flux front in the superconductor essentially depends 
on the anisotropy of pinning. In particular, in the case 
corresponding to the intrinsic pinning in high- Tc super- 
conductors, the front is a very sharp step, which should 
be taken into account in analyzing data of local magnetic 
measurements. We shall also show that under certain 
conditions two penetrating flux fronts can occur in an 
anisotropic superconductor. 

As usual, we consider here the cases when the charac- 
teristic magnetic field in the sample is sufficiently large 
such that the difference between the magnetic induction 
B and the field H may be disregarded. This condition 
is satisfied when jcd is much larger than the lower criti- 
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cal field Hd (otherwise, tiie so-called geometric barrieita 
must be taken into account). We shall thus express all 
the following equations in terms of the magnetic field H, 
related to the current density by the Maxwell equation 
j = V X H. 



II. MODEL AND ITS SOLUTION 

We consider an infinitely long strip of width 2ti; and 
thickness d, filling the space —w < x < w, —d/2 < z < 
d/2, i.e. we place the y axis of the coordinate system 
along the central line of the strip and the z axis along 
the external magnetic field Ha which is applied normal 
to the plane of the strip. The increasing applied field 
induces a sheet current J along y, which is related to the 
z component of the magnetic field in the plane z = by 
the Biot-Savart law, 



1 

2^ 



J{t) dt 



(1) 



Here and below all singular integrals are taken in the 
sense of the Cauchy principal value. The penetration of 
the magnetic flux into the superconducting strip is de- 
scribed by the following critical state equations: In the 
flux- free central region |a;| < b{Ha) one has 



(2) 



while in the region b{Ha) 
ready exists, one has 



< |a;| < w, where the flux al- 



\J{x)\^MH,{x)] 



(3) 



The position x = b{Ha) of the boundary separating the 
regions, is found by solving these equations. In Eq. (3) 
Jc{Hz) is the critical value of the sheet current- |_|At. 
present an exact solution of Eqs. (l)-(3) is knownac3'EJ 
only for the Bean critical state model where = const. 
Below we shall obtain the exact solution for the more gen- 
eral case when Jc{\Hz\) has the model form (see Fig. 1): 



Jc{Hz) = JcO 



for 0<Hz<H'^, 



for 



(4) 



Here 7 — [Jd — Jco)/H^] the three parameters Jd, Jco, 
and may have any positive value. 

As was mentioned above, in the case of thin supercon- 
ductors the dependence of the critical current density jc 
on the angle 9 between the local direction of the mag- 
netic induction and the normal to the strip plane can 
be taken into account if one considers this superconduc- 
tor as infinitely thin but with an Hz dependent sheet 
current. The model dependence described by Eqs. (4) 
corresponds to |the following ^-dependence of the critical 
current densityts shown in Fig. 1: 
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FIG. 1. Visualization of the dependence of the critical sheet 
current on the perpendicular magnetic field [ Jc(-ffz), Eq. (4), 
upper plot], equivalent to an out-of-plane anisotropy \jc{d), 
Eq. (5), lower plot]. The model has three independent posi- 
tive parameters, Jco, Jci, and all of same dimension. In 
this plot we put Jco = 1 and show two examples: Jci = 2 (in- 
trinsic pining, solid lines) and Jd = 0.5 (dashed lines), with 
Hi? = 0.3 (0.6) equivalent to 6*0 = arctan( Jco/2i^°) f» 60 (40) 
degrees. 

]c{e) = Jco/d for < < , 
JciO)^Jci/d for eo<9<n/2, (5) 

where tan = Jco/2H^. Thus, the case 7 > de- 
scribes intrinsic rPjinning by the CuO planes in high- Tc 
superconductorscj {jc peaks at 6* = n/2), whereas the 
case 7 < can be used to analyze pinning by columnar 
defects normal to the film {jc peaks at 6* = 0). In both 
these cases one can find two-dimensional solutions of the 
critical state equations for strips of small but finite thick- 
ness using the results obtained below and Eqs. (5,6,9-11) 
of Ref. |l|. 

Accounting for the symmetry of the sheet current, 
J{—x) = -~J{x), we seek the solution of Eqs. (l)-(4) 
in the form 



J{x) 



-[Jo{x) + Ji{x)] 



where 



Jo{x) = JcO , 
7 f \ 2JcO 



arctan 



(w^ — 6^) x^ 

11)2 Qji2 _ 



52 < < 
1/2 



x^ <b'^. 



(6) 

(7) 
(8) 



while Ji{x) is a new unknown function. The parameter 
b defines the position of the flux front, i.e., x = b is the 
point where Hz goes to zero. This parameter depends on 
Ha and must be determined together with Ji(a;). Both 
Jo{x) and Ji{x) (and the magnetic field below) are even 
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functions, which depend only on The function Jq{x) 
has the form of the exact solutionci to Eqs. (l)-(3) in the 
case when Jc — Jco and the external magnetic field is 
equal to 

Hb = Hcs arccosh(u' /b) 

where Hcs = Jco/t^- Using Eqs. (1), (6)-(8), the expres- 
sion for the magnetic field can be rewritten as 



'''Ji(yi) ds 



(9) 



where a is defined by the equality Hz{a) = H^, and 
Ho{x) is the siim of Ha and the field generated by the 
current Jo(a;),cl 



Ho{x) =Ha-Hf,.. 
Hoix) =Ha-Hb 



i/cs arctanh 



(x2-62)y;2ll/2 



0<x^ < b\ 



¥ <x^ < wi 



(10) 



(11) 



In Eq. (9) it was taken into account that Ji{x) dif- 
fers from zero only in the region < x^ < where 
H,{x) < Hi 

With the above formulas, the critical state equations 
take the following form: In the interval < < &^ one 
has 



Ho{x) 



1 A" Ji{\/s) ds 
2^ 



and in the region < x^ < a? we arrive at 

Mx) , 1 f<h{^)ds 



Ho{x)-Hl^ 



7 



2ti 



(12) 



(13) 



In deriving Eq. (13) we have expressed Hz{x) for b^ < 
x^ < a} in terms of .J\ (x) using the equality 



Hz(x)^Hl 



Ji{x) 



(14) 



that follows from formulas (3), (4), (6), (7). Eqs. (12), 
(13) are linear singular integral equations with Cauchy 
type kernel- The theory of such equations is well 
elaborated,E3 and hence we can find a, b, and Ji{x) for 
any given Ha- 

To do this, we introduce the following notations: 

1 7^1 
a = — arctan — , p = a , 

TT 2 2 



= a, 

— ^„2 



q;_ = a + 1, 

+2^a+ u2 ,2\0 



F±{t) = (a^-i^)"± \t^-b' 

and define the function f{t) by the equalities 

f{t) = -2Ho{t), 0<t<b, 
f{t) = 2 siuTTQ • - Ho{t)] , b<t<a, 



i.e., f{t) is discontinuous a,t t = b. Then, the solution 
of Eqs. (12), (13) can be represented as follows: In the 
interval < x^ < b^ one has 



w N 2, ,^ , , r f{t)dt 



while in the interval b^ < x^ < a? we arrive at 

7 



Ji {x) = cos Tra 



(15) 



(16) 



and Ji{x) — for a?' < x^ < . Here the integrals are 
taken in the sense of the Cauchy principal value; F+ and 
_F_ refer to positive and negative values of 7, respectively. 
If 7 < 0, for the above solution to exist it is necessary 
that 



"fit) 



dt^Q, 



and 



F^it) 



dt = 0. 



(17) 



(18) 



These two equalities enable us to determine b and a when 
7 < 0. If 7 > 0, the necessary condition for the existence 
of the solution is 



°/(^) 
oF+{t) 



dt = 0. 



(19) 



A second relation between a and b in this case is ob- 
tained from the analysis of the magnetic field near the 
point x'^ — a^. It turns out that 



Hzix) - Hi - C±^(4 + l'Y"{x- - a^r^ 
2|7l 

if tends to a? from above, and 

H,{x)-Hl^C±{a^-x^r^ 



(20) 



(21) 



if x'^ approaches from below. Here C± are certain in- 
tegrals independent of x] the subscripts -I- and — refer to 
the cases of positive and negative 7, respectively. Since 
Hz{x) > H° when x"^ > c? , we find that C+ > 0. On 
the other hand, one has Hzix) < H^ when < a^, and 
thus C± < 0. Hence, one concludes that C+ — 0. This is 
the second equality in the case of positive 7, and it has 
the form 



fit) dt 



/(«) 



fit) 



/(«) 



t(t2_62);3 a(a2„^2)/3 



tdt 



(a2-i2)i+« 



= 0. 



(22) 



The Eqs. (19) and (22) determine a and b when 7 > 0. 
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FIG. 2. Some profiles of the sheet current J{x) and of the 
perpendicular magnetic field Hz{x) in a superconductor thin 
strip with width 211; for various Jc{Hz) dependences, Eq. (4), 
equivalent to various out-of-plane anisotropies, Eq. (5), in an 
applied field Ha ~ 0.5. The unit for both J and H is Jco — 1. 
The anisotropy parameters are H° = 0.6 and Jci = 0.25, 
0.5, 0.75, 1, 1.5, 2.2, 4, and oo. The isotropic (or Bean) case 
Jci = 1 is shown as bold lines. The dotted lines indicate the 
field Hz = Hz and the position x = a, where J(a) = Jco and 
Hz{a) = _ff°. In the limit Jd — > oo the field Hz{x) at the 
flux front X — b abruptly jumps to the value Hz and stays 
constant for b < x < a. 



III. ANALYSIS 

Let us now analyze the obtained solution. For evalu- 
ation of the integrals in Eqs. (9), (15)-(19), (22) we use 
the method given in Appendix A. Some profiles J{x), 
Eq. (6), and Hz{x), Eq. (9), obtained in this way are 
shown in Figs. 2 to 5. 

It should be noted that no restriction on C_ is ob- 
tained when 7 < 0. In this situation the constant C_ is 
not equal to zero but negative, and thus the derivative 
of H-z with respect to x becomes infinite at x = a. In 
the same point a sharp bend occurs in J{x). In other 
words, we obtain that two flux fronts exist in the sample, 
at X = 6 and at a; = a, see Figs. 2, 4, and 5. Of course, 
the singularities in Hz and in J at x = a result from the 
sharp bend in our model Jc{Hz) at — H^, see Eqs. (4) 
and Fig. 1. However, one may expect that in the case 
7 < our qualitative conclusion on the existence of the 
second flux front in the sample remains valid if Jc{Hz) 
is a smooth function but its behavior changes abruptly 
over an interval smaller than Hcs- Such changes indeed 
may occur if the critical current density has sufficiently 
sharp angular dependence jc{d). 

We shall now describe Hz{x) and J{x) in the vicin- 
ity of the point x = b in which Hz ~ 0. According to 
Eq. (14), at this point \J{b)\ = J^. When x^ < 6^, it 
follows from the exact solution that 



\Jix)\-J,,^C^ib'-x^f. 
while if x'^ > b^, one has 



I J(a:)| - Jci 



(4- 



_7 

y)i/2 



ix' - b') 



(23) 



(24) 



Here and are certain integrals which do not de- 
pend on X and have negative values. Formulas (23) and 
(24) show that in the case 7 > 0, |J^(x)| has a sharp peak 
at X = &, whereas for 7 < 0, J{x) is a monotonic func- 
tion and its derivative with respect to x becomes infinite 
at X = 6, see Figs. 2, 4, and 5. Taking into account the 
above formulas and Eq. (14), one obtains the distribution 
of Hz near x ~ 6, 



Hz^Q 



for X <b , 



H, = -- 



(25) 
(26) 



,,-(x' -6')'^ for x>b. 
(4 + 72)1/2^ 

When 7 = 0, we arrive at the well-known result0@'0 
Hz oc (x^ — 6^)^/2. However, in the general case, tak- 
ing into account the equality = 5 ~ 7 arctan(7/2), one 
may conclude that the greater 7 is, the sharper is the 
Hz profile. Fig. 2. Interestingly, the dependence (x — b)^ 
sufficiently well describes Hz (x) even if x is not too close 
to 5, see Fig. 6. 

Consider now the solution in the limit of small positive 
values of 7. If 7 ^ 0, two cases are possible: Hz remains 
a constant, or it increases as 7'^ (i-e., Jci~ Jca — const.). 
In the first case one has a « 7/27r, Ha — iJf, oc 7, and the 
function f{t) tends to zero. Thus, according to Eqs. (15) 
and (16), Ji j=dJL_and the solution goes over to the well- 
known resultclElEj for the Bean critical state model with 
Jc — JcQ- In the second case Jo(x) -t- Ji(x) also tends 
to the solution corresponding to a constant Jc but now 

In the limiting case 7 +00, this parameter drops 
out from Eqs. (15), (16), (19), (22), and Ji(x) depends 
only on 7?°, Jco- In other words, if Jd » Jco,Hz or 
Jci > JcO ^ Hz J the solution becomes practically inde- 
pendent of Jci . The distribution of the magnetic field in 
this case can be understood using Eq. (26). It turns out 
jHz for 7 3> 1, and hence 



that 



Hz{x) « H°z{x^ - b 



(27) 



with /3 — !■ 0. This means we have an abrupt step of height 
H° at x^b (see Fig. 2). 

It should be emphasized that this limiting case, 7 — > 
-f-oo, corresponds to intrinsic pinning in high- Tc super- 
conductors in which the ratio [7c(7r/2)/7c(0)1 = Jci/Jco 
can be sufficiently large (see, e.g., Ref. ^J). Thus, our 
solution of this limit can be used for analyzing the crit- 
ical state in these superconductors. In particular, it 
follows from Eqs. (19), (22) that the position of the flux 
front, b/w, is a function of Ha/ Hcs and of the parameter 
Hz /Hcs, see Fig. 7. In general this function can not be 
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FIG. 3. Profiles of the sheet current J{x) (top) and of 
the magnetic field Hz{x) in a thin strip with width 2w and 
anisotropic pinning (solid lines) in an increasing applied field 
Ha = 0.15, 0.3, 0.5, 0.8, and 1.2 in units of Jco = 1- The 
anisotropy parameters are Jci/JcO = 1-5 and H^/JcO ~ 0.5, 
thus 7 = 1. The dashed lines show the profiles of an isotropic 
strip for the same values of the front position bo{Ha), Eq. (28). 
Note the sharp peak of J{x) at x — b of height J{b) = Jd and 
the steep front of Hz{x) at x = b for this type of anisotropy. 
At a: = a, J{x) reaches the value Jco = 1 and H^ix) goes 
through the value Hz{a) = H'^ marked by a dotted line. 

fitted by scaling the dependence found in the isotropic 
case,D 



6o 



1 



cosh{Ha/Hcs) 



(28) 



using some effective value of Hcs- Rather, the shape of 
b{Ha) essentially depends on the ratio H^/Hcs- There- 
fore, measuring b{Ha) in principle can give information 
not only on H^s — Jca/T^ but also on H^, i.e., about the 
width of the peak in jc{0), see Eq. (5). In particular, 
when H° < Jco, Eqs. (19), (22) lead to the following 
expression for the front position: 



1 + fc ■ (gP/g,,)'tanh^(gJg, 
cosh' {Ha/ H,,) 



(29) 



where the constant fc is determined by the root of the 
equation 
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FIG. 4. As Fig. 3, but for different type of anisotropy, 
Jci/Jco = 0.67 and ff^/Jco = 0.5, thus 7 = -0.67. In this 
case J(x) is monotonic and has an infiection point with ver- 
tical slope at X = 6 where J(6) = Jd (dotted line). The pen- 
etrating front of Hz{x) is now less steep than in the isotropic 
case, which is shown as dashed lines. 

— (u — 1) = u — arctanw , 
16 

fc= ^T. ^ « 0.394. 

7r2 (1 -I- 7/2)2 

Note that the right hand side of Eq. (29) cannot be re- 
duced to the dependence (28) in the whole interval of 
changes of Ha when i/^ is different from zero. The exact 
values of the front position b{Ha) are shown in Fig. 7 for 
the limit of large 7^1, for Jd = 11 and H^ = ... 1.5 
in units of Jcq — 1. 

In the third limiting case when "/ —00, one has 



/? ^ 1, 
comes 



-W^ and the the induction profile be- 



(30) 



with a small prefactor of the order of H^/\^\. Thus, for 
—7 ^ 1 the flux front at a; = 6 practically disappears 
while, according to Eq. (20), the second front near x — a 
is well developed, see Fig. 5. 
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FIG. 5. Profiles of J{x) (left) and H^{x) (right) in a thin 
strip with anisotropic pinning of the type Jci = 0.5 for various 
values of H° = 0.02, 0.2, 0.4 and 0.7 in a constant applied 
field Ha ~ 0.5 (in units of Jco = 1). The dotted linos at 
X = b, X = a, J = Jci, and = shall help to identify 
the features J(6) = Jci, J(o) = Jco, and Hz{a) = H°. Note 
that with decreasing the penetrating fiux front at x — b 
becomes less pronounced and a now front appears at x = a. 
In the limit only the front at a; = a remains and the 

profiles look like in the isotropic strip with b replaced by a. 

Finally, wc consider in some detail the case of small 
negative values of 7 when H° ^ Hcs = -/co/i" while the 
ratio Jco/Jci is not close to unity. This case can give 
some idea of pinning by columnar defects, which pro- 
duce a peak in jc{9) a^t 6 = 0. Indeed, if one assumes 
that the characteristic width of the peak, 9o, is small 
(^0 ^ 1)) then it follows from the definitions of H° and 
7 that « Jco/200 and I7I < 29o. Since the solution 
with 7 = and Jc = Jci describes the critical state in the 
strip before the irradiation [we assume that the colum- 
nar defects do not change jc{9) aX 6 > 9q\, the difference 
between the solutions corresponding to 7 7^ and 7 = 
provides information on pinning by columnar defects. In 
the considered case this difference is small, and it can be 
analyzed analytically. In particular, we obtain the fol- 
lowing relation between the positions of the flux fronts, 
h and hi, obtained at the same Ha in the strip with and 
without columnar defects, respectively: 
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FIG. 6. Comparison of the shape of the profile Hz{x) near 
the fiux front with the expression (a; — b)'^ suggested by 
Eq. (26). Shown are examples with Ha = 0.5 (in units of 
Jco ~ 1) and two different anisotropies; — 0.5, Jci = 1.5 
(thus 7 = 1, a = 0.148, P = 0.352, b = 0.477, a = 0.711) and 
H° = 0.5, Jci = 0.5 (thus 7 = -1, a = -0.148, p = 0.648, 
b = 0.213, a = 0.784). The exact H^{x) (dotted lines) is weh 
fitted over a large interval of x by the function c ■ (x — b)^ 
(solid linos) with c — 0.840 or c = 0.580 for these two exam- 
ples (with x and b in units of the strip half width w). The 
solid lines in the upper plot are straight lines fitting Hz(xY^^ . 

the exact dependence h{Ha) can not be described by 
Eq. (28) with some effective Hcs- The prefactor 



|7| ^ 2^0 je(0)-Jc(7r/2) 

in Eq. (31) is determined by the characteristics of pinning 
by the columnar defects, i.e., by the width and height of 
the peak in jc{9). 



(31) 



where h = nHa/Jci, w/bi = cosh(/i), and the function 
g{h) has the form: 



IV. CONCLUSIONS 



g{h) = [ ln{2cosht) dt. 
Jo 



Since (7 is a nonlinear function of h, 
g{h) « ^/i^ + 0.411 (1 - 



(32) 



(33) 



An exact solution of the critical state equations for 
the strip in perpendicular magnetic field is derived for 
an induction-dependent critical sheet current Jc{Hz) 
described by Eqs. (4). This model dependence may 
be used to simulate the intrinsic pinning by CuO planes 
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FIG. 7. The position h of the flux front, or penetration 
depth w — &, of a superconductor thin strip with width 2w 
and various Jc{Hz) dependences, Eq. (4), plotted versus the 
applied magnetic field Ha in units of Jco = 1- The dotted 
lines arc computed as described in Set. Ill for anisotropy 
parameters Jd = 11 and H" = 0, 0.3, 0.5, 0.7, 1, and 
1.5. The bold solid lines for = and 0.3 arc from 
Eq. (29) and fit the exact data very well. The dashed 
line h/w — l/cosh[Ji"a/(2.1i/cs)], obtained by stretching the 
isotropic (Jc = const) expression, Eq. (28), by a factor of 2.1, 
demonstrates that such scaling of the isotropic result cannot 
fit the anisotropic result. 

(7 > 0) or pinning by extended defects (7 < 0) in high- 
Tc superconductors. In the case 7 > 0, the pro- 
file in the vicinity of the flux front is sharper than in 
the isotropic case, and the current density has a sharp 
peak there. In the limiting case, 7^1, which may de- 
scribe the intrinsic pinning in high-Tc superconductors, 
the field profile Hz{x) has a sharp rectangular step. In 
the opposite situation, < 0, two flux fronts can oc- 
cur in the superconductor; the profile near x = b 
is less steep than in the isotropic case, and the current 
density is a monotonic function of x. In both cases 
of positive and negative 7 the profile Hz{x) in a suffi- 
ciently large vicinity of the flux front is well approximated 
by the expression Hz{x) w (x — b)^ with the exponent 
;3'= 0.5-7r-iarctan(7/2). 

The experimental investigation of flux-density profiles 
near the flux front and of the dependence of the pen- 
etration depth can give information on the strength and 
anisotropy of flux line pinning in superconductors. 
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APPENDIX A: NUMERICAL EVALUATION 

The condition that two integrals have to vanish, e.g. 
Eqs. (17,18) of the form /i(a,6) = and hia.b) = 0, 
we satisfy by minimizing the function U (a, b) = if + J| 
with respect to a and b. After this we calculate the sheet 
current J\{x) from Eqs. (15,16) and the magnetic field 
Hz{x) from Eqs. (9) and (14). 

The integrals (9), (15-19), and (22) over the variable 
t have integrands which possess one or several infinities 
at the points t = Q, t — x,t = b and t = a where the 
denominators vanish. We evaluate such integrals in the 
following way. 

In the integrals containing a factor {t — x)~^ we sub- 
tract the singular part and integrate it analytically, e.g., 

ri^dt=ri^^dt-M,^^_±ji. (Ai) 

Jo i — x'^ Jo ^ 2^ 2x a — x 

Then we divide the integration interval into pieces 
bounded by the remaining singularities, < t < b, 
b < t < a, and a < t < 1. In each interval we substi- 
tute the integration variable by an appropriate function 
t = t(v) and integrate over u such that the new inte- 
grand has no infinity and vanishes rapidly at the bound- 
aries. This new integral may thus be evaluated as a sum 
over an equidistant grid u, with constant weights. For 
example we write 

/ g{t)dt = / g[t{u)] t'{u) du « Y.9iWi (A2) 

with g, = g[t{ui)], u, = (?: - l/2)/N, w, = t'{tii)/N, 
t'{u) = dt/du, and i = 1, 2, 3, ... iV. This integra- 
tion method is very accurate if the substitution is chosen 
such that the weights and the products giWi vanish 
rapidly at the integration boundaries, e.g., Wi ^ and 
Wi ^ {\ — UiY with p S> 1 and g 3> 1. Simple choices of 
this substitution in the example (A2) arc 

t{u) = (3m2 - 2u^)t, t'{u) = 6m(1 - u)t , (A3) 

or better, 

t{u) = {10u^-15u^ + 6u^)t , t'(M) = 30m2(1-u)V. (A4) 

Higher accuracy is achieved by the following substitution. 
We chose equidistant Uj = (i — 1 /2)/N as above and then 
iterate (A3) m times starting with Si = Ui and Wi = t/N 
according to 

w := 6(s — s^)w , s := 3s^ — 25"^ (m times) . (A5) 

Finally we write t{ui) = SiT. The weights Wi = t'{ui)/N 
of this substitution vanish at the boundaries with expo- 
nents p = q = 2^'"^^), which can be made arbitrarily 
large. For example, using m = 5 iterations one gets the 
exponents p = q' = 2^ = 16. 

An infinity g{t) cx l/f^ in the original integral (A2) 
leads, after this substitution, to a new integrand vanish- 
ing at t = as g[t{u)]t'{u) oc with 1? = p{l — r]) — r]. 
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Thus, for the example 77 = 1/2 with p = 16 the new in- 
tegrand near u = vanishes as u"^'^ and the terms in the 
sum (A2) as {i — l/2)'^-^, in spite of the singular orig- 
inal integrand. For general exponent 7], to reach high 
accuracy one should choose m so large that the new ex- 
ponent is (1 — 77)2™^^ — 77 > 4, or approximately 
m > 3.5 — 1.5 ln(l — rj). To avoid spurious results due to 
rounding errors, one has to add in all vanishing denomi- 
nators a small e « 10^^^ by writing, e.g., Qt^ — b^\ + e)^ . 

In the limit of a large negative slope 7 ^ — 00 one 
has (3^1 and the integrals (17,18) containing a factor 
\t'^ — b'^\~^ are close to diverging. In this case the singular 
part in these integrals should be integrated analytically, 
similar as shown in Eq. (Al). The subtracted terms are 
conveniently chosen such that the integral which has to 
be taken analytically is simple, e.g., J t ■ {b^ — t^)~^dt. 
Note that the numerator f{t) in Eqs. (15-19,22) is dis- 
continuous at t = 6. 
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